CMB Bispectrum from Non-linear Effects during Recombination 
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We study cosmological perturbations by solving the governing Boltzmann and Einstein Field equa- 
tions up to second order, and calculate the corresponding CMB bispectrum during recombination. 
We include all the second-order Liouville and collision terms, truncating the multipole hierarchy 
at I = 10, consistently including all m 7^ terms when calculating the bispectrum in the flat-sky 
limit. At this stage, we focus on contributions at recombination and we neglect 2nd-order vector 
and tensor terms, lensing effects, and late time non-linear ISW. We find that the signal-to-noise for 
the bispectrum is 0.69 for Imax = 2000, yielding an overall signal _Fnl = 3.19 (normalised relative to 
the local model). We find that the effective /nl's of the equilateral and local type are 5.11 and 0.88 
respectively. This recombination signal will have to be taken into account in a quantitative analysis 
of the Planck data. 

PACS numbers: 98.80.-k, 98.80.Es 
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Introduction. With the Planck data shortly 
to become available, we may be able to detect 
non-Gaussianities in the Cosmic Microwave Back- 
ground(CMB) for the very first time. Such a detection 
could provide crucial insights advancing our understand- 
ing of the physics of the very early Universe, so its im- 
portance cannot be overstated. However, as well as pri- 
mordial effects when the initial conditions are laid down 
during inflation, there can be late-time non-linear inter- 
actions of the photon-baryon fluid with gravity(For more 
details, refer to fH). There have been several past at- 
tempts at estimatingthe CMB bispectrum from the lat- 
ter in the literature [2-1^1 • A detailed quantitative analy- 
sis of these late-time contributions is necessary, not only 
to debiase potential primordial signals, but also as an in- 
dependent test of GR and cosmological perturbation the- 
ory. In this paper, we present our numerical implemen- 
tation of the calculation of the bispectrum produced at 
recombination and beyond. We obtain /j^£''' — 0.88, in 
agreement with previous work focussing on the squeezed 
limit [13, 14 1 . We also find an equilateral signal with 
/atl*' =5.11 which agrees very well with The over- 
all signal-to-noise {S/N)i-cc = 0.69 provides an effective 
Fnl — 3.19 contribution which is larger than antici- 
pated previously and should be incorporated in the anal- 
ysis of the Planck CMB data with a forecast variance of 
Ai^NL — 5. First, however, we briefly describe the un- 
derlying analytic methodology before describing the nu- 
merical pipeline that has been developed, leading to the 
recombination bispectrum results. A much more detailed 
discussion of these methods and results will be presented 
in a longer paper shortly [l8| . 



Second- Order Perturbations and their Bispec- 
tra. To study non-Gaussianities generated by non-linear 



effects, we solve the 2nd-order Boltzmann Equations 

£[II] ^ £[I,I] ^ ^ (^[1,1] 

where £ is the Liouville operator, £ is the collision op- 
erator, the superscripts [II] and [1,1] denote linear terms 
of purely second-order perturbations and cross terms of 
two first-order perturbations respectively. In addition, 
we require knowledge of the first-order ionization frac- 
tion [6] at recombination, which encodes the effect of the 
perturbed last-scattering surface(LSS), when evaluating 
£[1,1]. 

We begin with the metric defined up to second-order 
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dxf'dx" = (1 + 2$)dT2 - 2.S§idx'dT 
-[(1 - 2*)% -t- 23^ij]dx^dx^ 



(2) 



where a(T) is the scale factor, and $, ^I', S§iS and ^ij 's 
are perturbations expanded in the following way 



(3) 



with A denoting the perturbation quantities above. We 
will use Newtonian gauge up to second-order through- 
out. Moreover, we ignore the Ist-order and 2nd-order 



vector and tensor perturbations, i.e. 







and 
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We compute the corresponding reduced bispectrum 
during recombination in the flat-sky and thin-shell ap- 



^ In general, the second-order vector and tensor perturbations are 
non-negligible as they are sourced by the first-order scalar per- 
turbations nonlinearly, however, we can expect these to be sub- 
dominant. 
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proximation 



' LSS 
(2^)2 



dkfdk^P{ki)P{k2) / dndr2dr: 



+1^3 + 2^3 (4) 

where LSS denoting the last-scattering surface, r = tq — t 
with the present conformal time tq, P{k) being the pri- 
mordial power spectrum, k^ being the component of k 
perpendicular to the tangent plane [l5| . The first-order 
source function S'l'l refers to s!^^ in while the second- 
order source function 52nd can be expanded into its lin- 
ear and quadratic parts S'2nd = S'^^^l + S*!^'^!. The lin- 
ear part 5*1^^1 refers to eq. (40) of 8] with the first-order 
perturbations replaced by the kernels of the correspond- 
ing second-order perturbationtQ. Analogous to well- 
known linear perturbation theory, S^^^^ contains contri- 
butions from intrinsic photon density, Sachs- Wolfe(SW), 
Doppler, integrated Sachs- Wolfe(ISW) effects of purely 
2nd-order perturbations which require the full set of so- 
lutions from the 2nd-order BoltzmannQ and Einstein 
Field equations flo|| (BE and EFE respectively). 

In contrast, the quadratic part S*!^'^! requires solutions 
of Ist-order perturbations only and can be read in con- 
figuration space a^ 

S'li'il = 2e-^[5*^9j-($ + *)|^ - ($ + ^)n'd,A 



4(-n^a,$ + l-OA + ($ - ^')n^a,$ 



-ful- 
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where S'^ is the screen projector 0, 4A = 51 /I is the 
fractional brightness, rt*'s are the components of the ob- 
servational direction, is the quadratic collision op- 
erator, and f is the optical thicknessQ 

We rewrite the 1st term in the 2nd line of Eq. ([5]) 
Art*9i$ — n'9i(A<i>) — $n'9iA and perform an integra- 
tion by part for the 1st term of r.h.s., to obtair0 



2e~nS'^d,(<P + *)^ - ($ + -iiWd^A 
o-n} 

4($' + ^')(A + $) ($ - *)7i'9,$ + 2*^''] 



2g[(4A + £)$ + 2$2 -I- C^^'^l] 



(6) 



where g = t' and £ = r'£. 



2 The kernels y^I"! (ki, k2, t) are defined as 
y^I"! (k, r) = f dkidka (27r)-3/25(k - ki - ka)^'"! (ki , ka, r). 

^ For simplicity, the notation for Ist-order perturbations is omit- 
ted. 

The quadratic collision CP'^l contains only quadratic terms in 
Ist-order perturbaions and refers to eq. (A. 55) of 4] in Fourier 
space. 

^ Another integration by part is needed to transform — n'di^)^ 
into ^{^' + (SW-ISW) and f'4>2 (SW-SW). 



In particular, the first and second term in the first line 
of eq. ([6]) are the lensing and time-delay effects respec- 
tively. The term A($' + ^f') can be identified as the 
photon-ISW coupling when the perturbed photons are 
redshifted by the time-varying gravitational potential. 
Similarly, we can interpret $($'-1-4''), (4A + £)$ and $^ 
as the SW-ISW, photon-SW and SW-SW effects respec- 
tively. The term ($ - \I')n'9j<i> -I- 24'^' is the quadratic 
part of the evolution equation of photon energy p in 2nd- 
order, i.e. {dp / dri)^^'^^ We will discuss the contributions 
of these effects on the bispectrum later. 

Our ultimate goal is to calculate the bispectrum as ob- 
served today, which means that all the effects post recom- 
bination till today are ideally included in the calculation. 
However, an appropriate milestone towards this goal is 
the computation of the bispectrum at the end of recom- 
bination which is what we undertake in this paper. This 
means that we cut off the numerical integration at this 
point. 

However, there exist source terms in this calculation 
which, via integrations by parts, are either zeroes or can 
be rescaled away today, i.e. we can drop the boundary 
terms if they are taken today. In moving from eq. ([5]) to 
eq. ([6]), as we mentioned, we have performed integration 
by parts for the term n^di{A^) in the line of sight integral 
and dropped the boundary term evaluated today. This 
means that we have secretly included post-recombination 
effects as a bonus. More importantly, these terms turn 
out to be much easier numerically to deal with once in- 
tegration by parts has been performed. 

Thus, it is necessary to clarify which terms of the 
source function are considered when we restrict the time 
integration around the recombination. In our study, we 
include all the terms of the source function S'2nd in eq. ([6]) 
coupled with the function g as it behaves like a dirac- 
delta function at recombination. We will also include all 
ISW-related effect at early-time0, such as the purely 2nd- 
order ISW and photon-ISW couplingQ. We ignore con- 
tributions from lensing, time-delay and late-time ISW- 
related effects. The len sing and time-delay effects have 
been studied elsewhere ll|, . Finally, we leave for a 
future publication the effects from the 2nd-order vector 
and tensor perturbations and the late-time ISW-related 
effects - note that the former contributes throughout the 
entire time integration of the line-of-sight approach al- 
though they should be fairly small. 

Numerical Implementation. To calculate the bis- 
pectrum in eq. (|3]), we solve the perturbations numeri- 



® We also take {dp/dr])^^'^^ into account but its contribution to the 
bispectrum is small. 

''' Effectively, we ignore the ISW effect due to late-time Dark En- 
ergy driven acceleration. We perform the time integration from 
T = 230 to 1050 Mpc to make sure that the early-time ISW 
effects are included. 
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cally for both 1st and 2nd-order terms. For the latter, 
we compute the kernels of the 2nd-order perturbations 
as functions of ki, ^2, ki-ki and r. 

The solved 2nd-order perturbations can then be fed 
into the purely 2nd-order source function 5'["1. Although 
only Ist-order perturbations are needed to compute the 
quadratic source function S'P'^', it has to be decomposed 
into multipoles sf'^j . In principle, to close the hierarchy 
requires multipoles up to infinite Z, but this is numeri- 
cally consuming, and hence the usual prescription is to 
truncate the hierarchy once convergence is reached. For 
the bispectrum at recombination, we found that I up to 
10 is sufficient (See Fig. [T]). We include all the m 7^ 
modes consistently up to / = 10. 




FIG. 1. The graphs of the bispectra generated from the 
quadratic source function S''''' against I for equilateral(upper) 
and squeezed(lower) limit. The curves correspond to differ- 
ent I truncations of the multipoles Sj^^' . We can see that the 
convergence occurs when I goes up to 10. 



chosen to align with the z-axis of the multipole decom- 
position. Thus, modes with m 7^ are suppressecJl. Sim- 
ilarly, the 2nd-order vector and tensor perturbations are 
not important in the squeezed limit as they are sourced 
by m — 1 and m ~ 2 modes respectively. 



Results and Discussion. In Fig. [2l we plot the bis- 
pectra of the terms in the source function S'2nd • We can 
see that the main contributions come from the effects of 
the Photon-SW, SW-SW and quadratic Collisions as well 
as the purely second-order SW, Doppler and anisotropic 
stresf@ effects. The sum of the Photon-SW and the SW- 
SW effects has roughly the same value regardless of I. 
That is, the total bispectrum from these two effects is 
approximately proportional to the product of the power 
of the long and short wavelength modes, i.e. Ci^Ci^. 
This is because the power spectrum of the short wave- 
length mode comes mainly from intrinsic intensity, SW 
and Doppler effects while the power spectrum of the long 
wavelength mode is proportional to the square of the ini- 
tial gravitational potential This shifts the bispectrum 
from the purely second-order source terms up and sup- 
presses the correlation between the bispectrum at recom- 
bination and that of the local type. 




Our bispectrum at the squeezed limit is in good agree- 
ment with the analytical estimate in 13|, ll^- However, 
as we have indicated previously, our integration by part 
scheme has included different post recombination effects 
and hence we do not expect a full agreement. 

In addition, we remark that the numerical accuracy of 
the multipoles sf'^ with m 7^ cannot be checked with 
the analytic solution as they contribute negligibly in the 
squeezed limit. The reason is as follows. In this limit, 
kj' <C k2 , k| , and combined with the conservation law 
of momenta, k2 is approximately parallel to k, which is 



* This is true if we only consider scalar perturbations in Ist-order. 
In this case, we have only m = modes in Ist-order when k 
aligns with the 2-axis of the multipole decomposition. 

^ Explicitly, the term 11 in . 
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FIG. 2. The bispectra of the terms in the quadratic(upper) 
and linear (lower) source function in squeezed limit. 

We present the full bispectrum at recombination in 
Fig. [21 showing isosurfaces of the bispectrum density. Al- 
though the main contribution to the recombination bis- 
pectrum is concentrated towards the edges of the tetrahe- 
dron, in the strongly squeezed limit the bispectrum fluc- 
tuates around zero and hence its correlation with the lo- 
cal type template is somewhat suppressed. What is clear 
from the figure is that the bispectrum does not corre- 
late particularly well with the popular templates - local, 
equilateral or orthogonal. It possesses its own distinct 
shape with most power located between the local and 
equilateral limits. In Fig. 21 we show different tetrahe- 
dral cross-sections through the full recombination bispec- 
trum taken at different summations h + h + h — const. 
On the other hand, it contains features in the squeezed 
limit and along the edges which reflect those appearing in 
the ISW-lensing bispectrum, but we will report on ISW 
cross-correlations in our longer paper [3l ■ 

In Table m we summarize the effective /nl's, the nor- 
malized Fnl's [i3 and their signal-to-noise ratios. The 
local-type F^£°-^ is small because of the fluctuating bis- 
pectrum in squeezed limit as we explained previously. 
The equilateral-type F^f^^' is small due to the lack of 
support in the interior of the bispectrum. 



TABLE I. The table of the effective /nl's, -Fnl's and signal- 
to-noise ratio S/N of the local and equilateral templates cor- 
related to the recombination bispectrum, as well as its total 
signal (auto-correlation). For ease of comparison, the _Fnl 
quantity normalises the integrated bispectrum signal for any 
shape relative to the /nl = 1 local model. It shows that 
the recombination bispectrum does not correlate well with 
local and equilateral templates. We have used Imax = 2000 
throughout. 



Model 


/nl 


Fnl 


(S/N) 


Equilateral 


5.11 


0.66 


0.028 


Local 


0.88 


0.88 


0.22 


Auto-correlation 




3.19 


0.69 



1200- 



800 




FIG. 3. The 3D plot of the reduced bispectrum generated 
at recombination. The bispectrum is normalized by the co- 
efficient D{h,l2,h) defined in [3l to remove an overall Z"* 
scaling. We can see that the main contribution occurs in a 
fairly diffuse region near the edges of the tetrahedron, i.e. 
towards the squeezed limit. The red regions represent pos- 
itive values while the blue regions represent negative values 
of the recombination bispectrum - the oscillatory nature of 
the bispectrum means that it does not conform well to known 
templates such as the local or equilateral types. The signal 
to noise for the bispectrum is 0.69 for imax = 2000 as plotted. 



These results are consistent in the squeezed limit with 
the correction of the value ~ 5 as 
. In this work, however, we have fo- 



14| . confirminj 
first calculated in 
cussed on increasing quantitative accuracy across the full 
tetrahedral domain beyond the squeezed limit. We have 
computed the equilateral type /nl"'' = 5-11 (-Fnl = 0.66) 
which is modest but consistent with However, the 
full signal-to-noise ratio of the recombination bispectrum 
in our calculation out to Imax = 2000 is substantial 
{S/N)rec ^ 0.69 (Fnl = 3.19) and larger than antici- 
pated in the recent literature. Confirmation of the mag- 
nitude of this signal will be important for the a quantita- 
tive analysis of the Planck data where this contribution 
should be incorporated in debiasing local and equilateral 
signals and in determining whether there is an overall pri- 
mordial non-Gaussian signature in the data. The recom- 
bination bispectrum will combine with the ISW-lensing 
bispectrum at about the 10% level (or greater) and its 
correlation can affect the significance of this determina- 
tion in the Planck data (recall that ISW can bias the local 
signal by as much as /nl = 9.5 [llj]). For this reason we 
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ism- 

FIG. 4. The cross sections of the normalized reduced bis- 
pectrum in Fig. |3]with the conditions + h + h) = 400 
(upper panel), 1100 (middle panel) and 1600 (lower panel). 
They correspond to hyperslices through Fig. |3] 



order (both which have been studied separately in the 
literature), the addition of the second-order vector and 
tensor perturbations, and finally the late-time ISW ef- 
fects. These will be addressed in a future publication [3l . 
With improving precision, this recombination bispectrum 
should be included in the analysis of future CMB exper- 
iments. 
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Note added in proof: While this paper was being pre- 
pared, a preprint appeared with an alternative numeri- 
cal implementation which also tackles second-order CMB 
perturbations [20| . This short paper offers a number of 
tests in the squeezed limit including a forecast /j^™' ~ 
0.82, in agreement with our results. An equilateral pre- 
diction was not given but the overall signal quoted 
{S/N)rec = 0.47 is lower than our result. This discrep- 
ancy may reflect the range of physical effects which have 
been incorporated in the two approaches. It is unlikely to 
be due to tensor and vector contributions which are sub- 
dominant, but it may reflect a different hierarchy trunca- 
tion scheme of the sf^ multipoles or the integration by 
parts approach they adopt which seems to include some 
different post recombination effects. 



shall continue to incorporate more physical effects in our 
numerical pipeline to improve this quantitative analysis 
further. 

Summary. In this paper, we focus on the bispectrum 
generated at recombination across the full range of multi- 
pole combinations. We find that the effective /nl's of the 
equilateral and local types are 5.11 and 0.88 respectively, 
while the overall signal-to-noise is {S/N)rcc — 0.69, all 
calculated using Imax — 2000. We note from Fig. [3] that 
the bispectrum possesses its own distinct features dif- 
ferentiating it from well-known templates, such as local, 
equilateral and ISW-lensing. To complete the full calcu- 
lation of this late-time bispectrum will require the inclu- 
sion of the time-delay and lensing effects up to second- 
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